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If quantum states exhibit small nonlinearities during time evolution, then 
quantum computers can be used to solve NP-complete problems in polynomial 
time. We provide algorithms that solve NP-complete and #P oracle problems 
by exploiting nonlinear quantum logic gates. It is argued that virtually any de- 
terministic nonlinear quantum theory will include such gates, and the method 
is explicitly demonstrated using the Weinberg model of nonlinear quantum me- 
chanics. 
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Computers are physical devices: like all other physical systems, their behav- 
ior is determined by physical laws. This seemingly obvious statement has im- 
portant implications, because it indicates that as our understanding of physical 
phenomena expands, the theoretical limits to the power of computing machines 
may grow accordingly. Recently, it has been shown that quantum computers can 
in theory exploit quantum phenomena to perform tasks that classical comput- 
ers apparently cannot, such as factoring large numbers in polynomial time|l|, 
searching databases of size M in time \fM Q , or simulating the detailed behav- 
ior of other quantum systems in less than exponential time and space || Q || . 
The realization that quantum mechanics could be used to build a fundamen- 
tally more powerful type of computing machine has led to a huge amount of 
recent activity in the field of quantum computation; for a review, see EkertQ 
or DiVincenzoQ. 

It has been suggested [§ (to) 0] (ill that under some circumstances the 
superposition principle of quantum mechanics might be violated - that is, that 
the time evolution of quantum systems might be (slightly) nonlinear. Such non- 
linearity is purely hypothetical: all known experiments confirm the linearity of 
quantum mechanics to a high degree of accuracy Q |lf| [jl6|] . Further, non- 
linear quantum theories have often been controversial and frequently have had 
theoretical difficulties |l7j [jl8]][|lj|. Nevertheless, the implications of nonlinear 
quantum mechanics on the theory of computation are profound. In particu- 
lar, we show that it is generally possible to exploit nonlinear time evolution 
so that the classes of problems NP and #P (including oracle problems) may 
be solved in polynomial time. An experimental question - that is, the exact 
linearity of quantum mechanics - could thereby determine the answer to what 
may have previously appeared to be a purely mathematical one. This paper 
therefore establishes a new link between physical law and the theoretical power 
of computing machines. 

The class NP is the set of problems for which, once given a potential an- 
swer, one can determine in polynomial time if the potential answer is in fact a 
solution. These include all problems in the class P (those that can be solved 
in polynomial time) as well as the NP-complete problems, e.g., the traveling 
salesman, satisfiability, and sub-graph isomorphism, for which no known poly- 
nomial time algorithms exist. We phrase our algorithm in terms of an oracle (or 
"black box"), which calculates a function that maps n bits into a single bit (i.e., 
it takes an input between and 2" — 1 and returns either or 1). We need to 
determine if there exists an input value x for which f(x) = 1; with a polynomial 
time algorithm to solve this problem, it is easy to solve all problems in the class 
NP. Indeed, this oracle problem is in fact a harder problem than those in NP, 
because it clearly requires exponential time on a classical Turing machine. 

A simple algorithm that solves the NP oracle problem can be thought of as 
an extension of Grover's data-base search algorithm^) to a nonlinear regime. 
Suppose that it is possible to perform a nonlinear operation on a single qubit that 
has a positive Lyapunov exponent over some (not exponentially small) region: 
i.e., somewhere on the unit sphere there exists a line of finite extent along which 
application of the operation causes nearby points to move apart exponentially 
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at a rate e A proportional to their original separation AO. We can exploit this 
behavior to solve NP problems if we begin with an ordinary quantum computer 
(i.e., one that can perform the usual quantum logic operations such as controlled 
rotation gates) and use the algorithm described below: 

Step 1. Begin by performing a it/ 2 rotation on each of the first n qubits to 
obtain the state 

Step 2. Use the oracle (only once) to calculate f(i): 



^iEi'./w) ( 2 ) 
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Step 3. Reverse the ir/2 rotation on each of the first n qubits. Each state 
|i>then maps into a superposition over all possible |i>, with amplitude ±"^pr- In 
particular, each state |i>contributes + -^=of its amplitude to the state |00...0>, 
for a total contribution of amplitude from each |i>. At least ^2™ of these 
states correspond to a particular value of f(i)=a, and thus the state |00....0,a>has 
amplitude at least 1/2. A measurement on the first n qubits will therefore yield 
the state |00...0>with probability at least 1/4. The system will then be in the 
state 
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where s is the number of solutions i for which f(i)=l. The last qubit now 
contains all the information that we need; however, for small s, a measurement 
of the last qubit will almost always return |0>, yielding no information. We 
wish to distinguish between the cases s=0 and s>0. 

Step 4. Repeatedly apply the nonlinear operation to drive the states rep- 
resenting these two cases apart at an exponential rate: eventually, at a time 
determined by a polynomial function of the number of qubits n, the number of 
solutions s, and the rate of spreading (Lyapunov exponent) A, the two cases will 
become macroscopically distinguishable. 

Step 5. Make a measurement on the last qubit to determine the solution. 
If the angular extent of the nonlinear region is small, it may be necessary to 
repeat the algorithm several times in order to determine the solution with high 
probability. In general, the algorithm will require 0((n/rj) 2 ) trials, where i] is 
the angular extent of the nonlinear region. If 77 is sufficiently large, the oracle 
may need to be called only once. 

To solve problems in the class #P, we need to determine the exact number of 
solutions s. This is approximately found by counting the number of times that 
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the nonlinear operator was applied. To determine s exactly, one proceeds with 
finer and finer estimates by rotating the final qubit such that the current best 
estimate is centered in the nonlinear region; in this way applying the nonlinear 
operator separates states with s near this value so that they are distinguishable. 
With only a polynomial number of iterations, one determines the value s ex- 
actly. The previous algorithm can therefore by used to solve arbitrary problems 
(including oracle probems) in the class #P as well. 

The above algorithm has one disadvantage in that it requires exponential 
precision. We describe below another algorithm, which (at least in the case 
of the Weinberg model) appears to be robust against small errors. In order to 
simplify the description, we assume (for now) that there is at most a single value 
x for which f(x) = 1. We begin as before with a quantum computer that can 
perform the usual quantum logic operations, and that can in addition perform a 
simple nonlinear operator whose form will be described below. We conceptually 
divide the qubits into n input qubits (representing the input, an integer between 
and 2™ — 1 ) and one flag qubit. In order to solve the stated oracle problem, 
we use the following algorithm: 

Step 1. We begin by performing a n/2 rotation on each of the first n qubits 
to obtain the state 

^iEM) ( 4 ) 
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Step 2. Use the first N qubits as the input to the oracle and store the output 
in the (N+l) s * qubit (the flag qubit). (Note that this is the only time that the 
algorithm needs to call the oracle; hence, we require only 1 evaluation of the 
function.) The system is now in the state 

^iE^w) ( 5 ) 

Step 3. Consider the first qubit separately, and group all the states of the 
superposition into pairs based on the value of qubits 2..n. That is, the qubits 
2..n define 2™ _1 subspaces of dimension 4 = (2 dimensions for qubit 1) * (2 
dimensions for the flag qubit). Within each subspace, the computer will be in 
one of the following states (where we write the value of the first qubit followed 
by the value of the flag qubit, and ignore the normalization constants): 

(a) |00) + |11) 

(b) |01) + |10) (6) 

(c) |00) + |10) 

(At the start of the computation, most of the superposition will be in the 
third state, because the flag qubit is |l>in at most only 1 of the 2™ components.) 
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A distinctly nonlinear transformation "N" is then applied to these two qubits 
(we show below how virtually any deterministic nonlinear operator can be recast 
into this form): 



(a) |00} + |11) — |01) + |11> 

(b) |01) + |10) — > |01) + |11) (7) 

(c) |00) + 1 10) — ► |00) + 1 10) 

This transformation is like an AND gate - it ignores the index qubit and 
places the flag qubit in the state |l>if and only if either of the original compo- 
nents had the state |l>for the flag qubitj^j 

Step 4. The previous step is then repeated using each of the first n qubits 
as the index (and the remaining n-1 qubits to define the 2™" 1 subspaces). After 
each iteration, the number of components in the superposition that have a |l>for 
the flag qubit doubles. After n iterations, the flag qubit is no longer entangled 
with the first n qubits: it is either in the state |l>for every component of the 
superposition or the state |0>for every component of the superposition. 

Step 5. Measure the flag qubit to determine the solution. 

Thus, if one can perform the two qubit nonlinear transformation N one can 
find the answer to an NP-complete problem with certainty in polynomial (in 
fact linear) time, and using only a single evaluation of the oracle. It may be 
objected that the nonlinear operator N appears arbitrary and unnatural: indeed, 
it was selected exactly so as to be able to solve the stated problem. However, the 
apparently arbitrary operation N can be built using ordinary unitary operations 
and much simpler and more 'natural' single qubit nonlinear operators (that is, to 
the extent that any nonlinear operation in quantum mechanics can be considered 
'natural'). One possible technique for generating the transformation would be 
to use the following steps: first, act on the two qubits with the unitary operator 
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This transforms the states above as follows 
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(8) 



2 There is one subtlety regarding nonlinear quantum mechanics which we should address 
here. When the superposition principle is abandoned, it is not immediately clear how entan- 
gled qubits will evolve. We follow the Weinberg model, in which the time evolution for a joint 
system composed of two subsystems is specified in terms of a preferred basis of vectors for 
the tensor product Hilbcrt space. For the purpose of using the nonlinear dynamics to perform 
quantum logic, we specify the joint dynamics in terms of the basis {|b>}={|0...00>, |0...01>, 
. . . , |1..11>} for each subsystem. The Weinberg prescription is as follows: write the joint 

state for the system |^ r ) 12 as Yl a b l 1 /'i>)2 an d then act on each |i/>i))2 independently with 

b 

the nonlinear transformation N. 
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(a) -L [|00) + |11)] — |00) 

(6)-L[|01) + |10)]^|01) (9) 

(c) ± [|00) + |10>] — 1 [|00) + |01) - |10) + |11)] 

Next, operate on the second qubit with a simple one qubit nonlinear gate 
ra_ that maps both |0>and |l>to the state |0>. Thus 

(a) -L [|00) + |11)] |00) 

(6)_L[|01) + |10)]^|00) (10) 

(c) -L[|00> + |10>] — |A) 

The third final state is unknown because we have not bothered to specify 
how the non-linear gate acts on the state |00>+ |01>- |10>+ |11>. This omis- 
sion thereby allows for flexibility in choosing the gate n_. Whatever the state 
|A>may be, we can perform a unitary operation that will transform the first 
qubit into the pure state |0>whilc leaving the state |00>in place. The computer 
is then in one of the following states 



(a) |0)|0) 

(b) |0)|0) (11) 

(c) |0)(z|0) + y\l)) 

A second non-linear gate n + is now required that will map the state x|0>+ 
y|l>to the state |l>(for the particular values of x and y which result from the 
above steps but not necessarily for arbitrary x and y), while leaving the state 
|0>unchanged. After this gate is applied, the transformation resulting from the 
steps described so far is then: 

(a) -^[|00> + |11>] — ►|00) 

(6)-L[|01) + |10)]^|00) (12) 
(c)-L[|00) + |10)]^|01) 

The two qubit transformation N is then easily obtained with a NOT gate on 
the second qubit and a tt/2 rotation on the first qubit. 
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Having thus shown how to generate N, the question is now reduced to that of 
generating the simpler single qubit gates n_ and n + . If one considers the state 
of a qubit as a point on the unit sphere, then all unitary operations correspond 
to rotations of the sphere; and while such rotations can place two state vectors 
in any particular position on the sphere, they can never change the angle be- 
tween two state vectors. A nonlinear transformation corresponds to a stretching 
of the sphere, which will in general modify this angle. The desired gates n_ and 
n + are two particular examples of such operations. Excepting perhaps certain 
pathological cases (e.g., discontinuous transformations), it is evident that virtu- 
ally any nonlinear operator, when used repeatedly in combination with ordinary 
unitary transformations (which can be used to place the two state vectors in an 
arbitrary position on the sphere), can be used to arbitrarily increase or decrease 
the angle between two states, as needed to generate the gates fi_ and n + . Wc 
describe in detail how these gates can be obtained using the model of nonlinear 
quantum mechanics put forth by Weinberg. 

In Weinberg's model, the "Hamiltonian" is a real homogeneous non-bilinear 
function h(ip,ip*) of degree one, that isQ 

i'kTr— = ^kjrn; = h (13) 
and state vectors time-evolve according to the equation 

at &r k K ' 

Following Weinberg |J , one can always perform a canonical homogeneous trans- 
formation such that a two-state system (i.e., a qubit) can be described by a 
Hamiltonian function 

h = nh(o) (15) 

where 

n= |Vi| 2 + kH 2 (16) 



\H 2 



(17) 



It is easy to verify his solution to the time dependent nonlinear Schrodinger 
equation dll), which is 

Mt) = c k e" l ^ t (18) 

where 
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Wi(a) = h(a) — ah! {a) 
w 2 (a) = h(a) + (1 — a)h'(a) 



(19) 
(20) 



For nonlinear h(a), one sees that the frequencies depend on the magnitude 
of the initial amplitude in each basis state. Intuitively, one can imagine a 
transformation on the unit sphere which, instead of rotating the sphere at a 
particular rate, twists the sphere in such a way so that each point rotates at a 
rate which depends upon its angle 9 from the axis (clearly, this transformation 
involves stretching of the surface). One can exploit this stretching of the sphere 
to build the gate n_ as follows: 

Step 1. Perform a rotation on the first qubit by an angle <f> < 45°: 

|0) — >cos(0)|O) -sin(0)|l) (21) 
|1) — >sin(0)|O)+cos(0)|l) (22) 

Step 2. Time-evolve the system according to the nonlinear Hamiltonian 
h = nh(a). Thus 

|0> — ► cos(0)|O) - sin(0)|l> — > acos(0)|O) - /3sin(0)|l) (23) 
|1) — > sin(0)|O> + cos(<£)|l) — >jcos(4>)\0) +6sm((f>)\l) (24) 

where a, /3, 7 and 5 are phase factors. Because the initial amplitudes of the 
basis states are different in the two cases, the nonlinear Hamiltonian will cause 
the components to evolve at different frequencies. As long as these frequencies 
are incommensurate, there is a time t at which a=^=5=l and (3=-l (to within 
an accuracy e). (Further, this time t is a polynomial function of the desired 
accuracy e.) The net result of these two steps is then 



|0) — >cos(0)|O)+sin(0)|l) (25) 

|1) — > sin(^)|0) + cos(0)|l) (26) 

Step 3. Reverse the first step. Thus 

|0) — ► cos(20)|O) + sin(20)|l) (27) 

|1) — |1> (28) 



Essentially, we have reduced the angle between the two states by an amount 
2(f>. By suitable repetition of this procedure (that is, by choosing (f> appropriately 
for each iteration), or simply by choosing <j> precisely in the first step, the states 
|0>and |l>can be mapped to within e of the state |0>, in an amount of time 
which is a polynomial function of the desired accuracy. This is the desired 
behavior for the nonlinear gate n_. The procedure can be modified slightly 
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to increase the angle between state vectors and produce the desired behavior 
for the gate n + . We have thus shown explicitly how to solve NP-complete 
problems using the Weinberg model, using an algorithm which did not require 
exponentially precise operations. 

To solve the problems in the class #P, one replaces the flag qubit with a 
string of log^n qubits and modifies the algorithm slightly - so that it adds the 
number of solutions in each iteration rather than performing what is effectively 
a one bit AND. In this case, a measurement of the final result reveals the exact 
number of solutions. 

We have demonstrated that nonlinear time evolution can in fact be exploited 
to allow a quantum computer to solve NP-complete and #P problems in poly- 
nomial time. We have shown explicitly how to accomplish this exponential 
speed-up using the Weinberg model of nonlinear quantum mechanics. In con- 
cluding, we would like to note that we believe that quantum mechanics is in 
all likelihood exactly linear, and that the above conclusions might be viewed 
most profitably as further evidence that this is indeed the case. Nevertheless, 
the theoretical implications and practical applications that would result from a 
discovery to the contrary may warrant further investigation into the matter. 
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